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Abstract 

We evaluate the complete four-loop anomalous dimension matrix that is necessary for 
determining the effective flavour-changing neutral current couplings qq >r y and qq'g at the next- 
to-next-to-leading order in QCD. The resulting 0{a 2 s {^)) correction to the B — > X s 7 branch- 
ing ratio amounts to around —2.9% for \i h = 5 GeV, and —4.4% for fi b = 2.5 GeV. 



1 Introduction 



Flavour- Changing Neutral Current (FCNC) decays are well known to provide stringent con- 
straints on extensions of the Standard Model (SM) because their SM amplitudes are loop- 
suppressed. An additional chirality suppression factor m q /M w for radiative FCNC transitions 
q — > q >r y makes these processes particularly attractive. The chirality suppression may be off- 
set in certain new-physics models like the Minimal Supersymmetric Standard Model with large 
tan j3 [1] or left-right models [2]. Precise SM calculations of the radiative flavour-changing decay 
rates are thus vital for our ability to derive constraints on new physics from these observables. 

In the present paper, we evaluate the complete four-loop Anomalous Dimension Matrix 
(ADM) that is necessary for calculating the effective FCNC couplings qq'^y and qq'g at the 
Next-to-Next-to-Leading Order (NNLO) in QCD. Our main motivation is contributing to the 
NNLO calculation of the inclusive B — > X s 7 branching ratio that is generated at the quark 
level by the b — > 57 transition. For this decay, both the current experimental errors and the 
non-perturbative effects are smaller than the perturbative NNLO corrections [3]. Our four-loop 
0(al(nb)) contributions turn out to suppress B(B — > X s ^) by around 2.9% for \i h = 5 GeV, 
and 4.4% for //& = 2.5 GeV. The dominant part of this effect has already been included in the 
phenomenological NNLO analysis of Refs. [3,4]. 

The article is organized as follows. In Section [2J we introduce the relevant effective La- 
grangian and describe the structure of the Renormalization Group Equations (RGEs) for the 
Wilson coefficients. In Section [31 the bare four-loop calculation is described. Section 0] is de- 
voted to the calculation of subdivergences. Our final results for the renormalization constants 
and ADMs are given in Section [51 The numerical size of the four-loop effects is analyzed in 
Section El We summarize in Section [71 Appendix A contains solutions to the RGEs for all 
the Wilson coefficients that matter for b — > 57 at the NNLO before including higher-order 
electroweak corrections. 



2 Effective Lagrangian 

All the FCNC processes that have been observed so far are generated at the electroweak scale 
fiQ ~ Miy,m t , while the relevant momentum transfer scale /i/ is much lower. Large logarithm 
{a s hi jjL 2 ^j resummation at each order of the perturbation series in a 8 (jif) is necessary to 
obtain viable results. It is achieved by making use of the RGEs in an effective theory that 
arises from the SM after decoupling of the heavy electroweak bosons and the top quark [5]. 

The effective Lagrangian that matters for the q — > g'7 and q — > q'g transitions at the 
leading order in the electroweak interactions is given below. For definiteness, we specify the 
quark flavours as in the b — > 37 case. 

4Gf 8 

£eff = ^QCDxQED^, d,8,c,b) + -^-V'Vt^CiiriQi + .... (2.1) 

V2 i=1 

Here, Gp and Vij stand for the Fermi coupling constant and the quark mixing matrix elements, 
respectively. The Wilson coefficients Ci(fi) play the role of coupling constants at the effective 
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weak vertices (operators) Qi. The operator basis is chosen as in Ref. [6] to avoid difficulties 
with 75 in dimensional regularization. 

Q x = (s lt ,T a P L c)(crT a P L b), 

Q2 = (s llx p L c)(crp L b), 

Q 3 = (s llx P L b)E q (qrq), 

Q 4 = (s^T a P L b)E q (qrT a q), 

Q 5 = {s lflll ^^P L b)Y, q {qi^Y"q), 

Q 6 = (s^^ 3 T a P L b)j: q (q^^T a q), 

Qr = jjfct[s<r'°'(rrisPL + m b P R )b]F lu ,, 

Qs = ^[sa^{m s P L + m b P R )T%]G% 

where Pl,r = (1 =F 7s)/2. For simplicity, terms proportional to the small V ub mixing have been 
neglected here. The dots on the r.h.s. of Eq. (12. ip stand for ultraviolet (UV) counterterms. 
Apart from Q\,. . . ,Qs, the counterterms contain unphysical operators that matter only off-shell 
and/or in D 7^ 4 spacetime dimensions (see Section H]). 

Instead of the original Ci(fi), it is more convenient to work with certain linear combinations 
of them, the so-called "effective" Wilson coefficients 

{Ci(fx), for i = 1, ...,6, 

CM/i) + £, 6 =1 ^(/i), for * = 7, (2.3) 
CM + Y.Uz 3 C^), for * = 8. 

The numbers yj and Zj are defined so that the b — > 57 and b — > sg amplitudes at the Leading 
Order (LO) in QCD are proportional to C7 and Cg ff , respectively [7]. In the MS scheme with 
fully anticommuting 75 which is used in our calculation, one finds, y = (0, 0, — |, — |, — y , — 
and z = (0,0,1, -|, 20, -f) [6]. 

The RGEs for Cf(/i) 

M^Cf (//) = CfMigM (2.4) 
are governed by the ADM 

fftfa) = ^ 5 ™+l 7 (n)eff (2.5) 
n>0 

where a = a s (n)/(47r). This matrix is determined from the effective theory renormalization 
constants. The matrices ^ n ^ eS have a block-triangular structure 



,(n) R (n) 
02x6 C*2x2 



7 (n)eff = ( ? x6 X) I • ( 2 - 6 ) 



The down-left block vanishes because the dipole operators Q7 and Q$ are actually dimension- 
five ones, and thus generate no UV divergences in dimension-six four-quark amplitudes. The 
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Figure 1: One of the O(10 4 ) calculated four-loop diagrams. 



diagonal blocks and C^ n > are found from [n + l)-loop renormalization constants. Non- 
vanishing contributions to the off-diagonal blocks B^ arise at two and more loops only. Once 
the normalization conventions for Q 7 and Q% are chosen as in Eq. (EQ]) . the blocks B^ contain 
information on [n + 2)-loop renormalization^] 

The complete matrices <y(°) eff and ^ 1 > eS together with the relevant references can be found 
in Ref. [6]. We quote these results in Section [5j The three-loop block B^ l > was confirmed 
in Ref. [8]. At the NNLO, one needs to know the full 8x8 matrix ^ 2 "> eS . The three-loop 
blocks and were calculated in Refs. [9] and [10], respectively. In the present paper, 
we evaluate the four-loop block B^ 2 \ 



3 Bare Four-Loop Calculation 

The renormalization constants contributing to the matrix B^ are found after sub loop renor- 
malization from the UV-divergent parts of four-loop diagrams like the one shown in Fig. [IJ 
Each of the operators Qi, ...Q% must be considered in the effective four-quark vertex, and 
the external gauge boson can be either a photon or a gluon. The overall number of relevant 
four- loop diagrams turns out to be 21986, when different color and gamma matrix structures 
are treated together. 

In order to simplify the calculation as much as possible, an infrared (IR) rearrangement as in 
Refs. [11-14] has been applied. In this approach, a Taylor expansion in the external momenta 
is performed after introducing a common mass in all the propagator denominators. The order 
of the expansion is determined by the dimensions of the operator insertion and the Green's 
function. In the present case, up to two derivatives have to be applied. Afterwards, the 
problem is reduced to evaluation of single-scale fully massive vacuum integrals. 

1 Inverse powers of coupling constants are sometimes used in the dipole operator normalization to make 
j( n ) cS a purely (n + l)-loop object. Such a choice is convenient at intermediate steps (see Section [5]), but the 
final results for the RGEs are more compact when the normalization as in Eq. (|2.2p is applied. 
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In principle, the contribution of any graph can be mapped onto the following Dirac structures 



where p and k stand for the incoming 6-quark and the outgoing photon/gluon momenta, re- 
spectively, whereas M denotes the regulator mass. However, due to the huge computational 
resource requirements, only the coefficients at the structures Sj for j = 7,9, 10, 11 have been 
calculated. Only these four structures are needed to determine the sought ADM (see 
Section SJ). 

Comparing the vertex diagram structure and the Taylor expansion depth with those needed 
in the evaluation of the four-loop /3-function with a single gauge parameter, we note that the 
maximal powers of numerators and denominators in the integrals are the same. This has allowed 
us to use the solution of the integration-by-parts identities and the master integrals determined 
in Ref. [14]. All the salient details of the employed techniques can be found in that paper. 

4 Subdivergences 

Considering renormalization of all the operators up to dimension six that can arise in the 
effective theory (12.11) is necessary for proper subtraction of subdivergences in our calculation. 
Since the four-loop diagrams are evaluated in the usual Feynman-'t Hooft gauge (without 
background fields), in D = 4 — 2e dimensions and off-shell, the following types of counterterms 
must be taken into account apart from the usual QCD ones: 

• the original ("physical") operators from Eq. ( 12.21) . 

• gauge-invariant operators that vanish by the QCDxQED Equations Of Motion (EOM), 

• gauge-variant EOM-vanishing operators, 

• BRST-exact operators, i.e. operators that can be written as Becchi-Rouet-Stora-Tyutin 
variations of other operators, and 

• evanescent operators, i.e. operators that vanish in D = 4 spacetime dimensions by Dirac- 
algebra identities, but have to be included for D^4. 

One begins with writing down all the possible operators of dimension up to six with proper 
flavour quantum numbers that fall into the above five classes. The potentially infinite set of 
evanescent objects is reduced to a finite one at each given order of the perturbation series once 
the prescription of Ref. [15] is applied. 

The long list of relevant operators can be further reduced by taking into account that the SM 
Lagrangian is invariant under the following discrete transformations of its fields and parameters 
(already after the electroweak symmetry breaking and mass matrix diagonalization): 

• b<->s, m b <-> m s , V Qb ^V Qs (Q = u,c,t). 




(3.1) 
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• CP transformation and V — > V*. 

A combination of these two transformations maps the AB = —AS = 1 interaction terms onto 
themselves, so it must leave C e s (12. ip invariant. The relative plus sign between the m s and m b 
terms in the dipole operators Qj and Qg is fixed by this very symmetry. Once the operators 
are determined, terms proportional to m s are often neglected. Their effect on the B — > X s ~f 
branching ratio is negligible (C(mJ/m^)) but they matter for the CP asymmetry. Some details 
on selecting the right number of gauge-invariant operators have been given in Section 5.3 of 
Ref. [16]. 

Explicit expressions for all the operators that satisfy the above-mentioned conditions can be 
found in Refs. [8,9], so we do not repeat them here. Let us only note that Section 2 of Ref. [8] 
contains all the EOM-vanishing operators, including the ones suppressed by the QED coupling 
e, and the one that gets generated at three loops only (Q22)- On the other hand, the only 
possible BRST-exact operator and all the evanescent operators that we need here are given 
in Section 3 of Ref. [9]. The BRST-exact operator remains irrelevant. In particular, we have 
checked that the three-loop analogue of the two-loop amplitude described in Fig. 3 of Ref. [9] 
turns out to mysteriously vanish, as well. 

Apart from all the "real" off-shell counterterms that have been mentioned so far, we en- 
counter certain spurious ones that arise due to our particular IR rearrangement. The dimension- 
three gluon-mass-like counterterm was discussed in detail in Refs. [11,12]. Spurious dimension- 
four flavour-changing operators that are proportional to the IR-regulator mass (see Eq. (18) of 
Ref. [8]) are included, too. 

The completeness of the full operators basis has been verified in many ways. For instance, 
since the number of operators contributing at the tree level to the b — > 57 transition is smaller 
than the number of independent Dirac structures Sj (13.11) . there exist seven linear relations 
between coefficients at these structures in the renormalized amplitude. All these relations 
together with three analogous ones for b — > sg have been successfully tested up to three loops [6, 
8,9,17]. In the present calculation, we determine coefficients at four of the structures Sj only, 
which leaves us with just one linear relation for b — > S7 and none for b — ► sg. Therefore, we 
essentially rely on the fact that the operator basis is indeed complete. 

Once the (n < 3)-loop diagrams with counterterms are added to the bare four-loop ones, we 
obtain an expression of the form 

Arfpp + m b (Agft^ + A 10l ^ + A U ^7 M ) (4.1) 

where the coefficients Aj contain UV-poles only. These poles need to be canceled by tree-level 
counterterms. First, we subtract such tree- level counterterms that originate from products of 
the usual QCD renormalization constants with the (n < 3)-loop operator ones, which modifies 
the coefficients in Eq. (14. ip . Let us call the new coefficients A\. The linear combination 

-A' - -A' + -A' - -A' (4.2) 

4 7 2 9 -r 2 10 2 11 v ) 

gives us the sought four-loop contributions to the renormalization constants and (for 
% = 1, ... 6) in the photonic and gluonic cases, respectively. The single consistency relation that 
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is successfully tested in the photonic case reads A' g + A' 10 + A' n = 0. Proving that the projection 
( 14.21) is the right one requires studying Feynman rules for each operator in the full basis and 
verifying that only Qj or Q$ can contribute to this particular combination at the tree level. 

A very powerful test of the results in Eq. (14.1 ft is provided by the so-called locality constraints. 
Let us denote the n-loop l/e m -contribution to Ai by v4 4 - nm ' ) , i.e. 

4 2ne 

Ai= £ ^r4 (nm) - (4.3) 

n,m=l 

The quantities A [ nm ^ must satisfy many relations to ensure that all the (ln fe //) / e J poles cancel 
out in Aj. Such a cancellation is necessary for the counterterms to remain polynomial in external 
momenta, i.e. local in the position space. The explicit relations read 



124 (44) = -34 (14) = 2A? A) = -34 (34) (4.4) 

6A (43) = - A (23) - 3 A (33) (4.5) 

g A (43) = A (13)_3 4 (33) (4 6) 

4A (42) = -A (12) -24 22) -3A (32) - (4.7) 



All these relations have been checked to hold in the photonic and gluonic cases for each of the 
six operator insertions. One should realize that testing the 1/e 2 poles in a four- loop amplitude 
constitutes a much more sensitive cross-check of the calculation than doing the same in a two- 
loop amplitude. Since Af m and Aj m have been evaluated with the help of different routines 
written by different authors, we can be practically sure that our calculation is free of errors like 
missing a diagram or taking it with a wrong overall factor. 



5 Renormalization Constants and ADMs 

The counterterms in the effective Lagrangian that we have calculated can be written in the 
following form 

AC 6 

ACtr* = -£-V*V tb £ Q (Z l7 X 7 + Z l8 X s ) , (5.1) 
V2 i=1 

where 

X 7 = 7^^Z Q 7 and X s = -^^Z Z g Z 1 G 2 Qs- (5-2) 

K z g9) ( Z a9) 2 

Here, Z , Z g and Zq are the usual QCD renormalization constants of the mass term, gauge 
coupling and the gluon field 

(m^) B = Z Q mW, g B = Z g g, G a / = Z^G^. (5.3) 

Below, we give our results for Z^ and Zjg up to four loops. The "trivial" terms proportional 
to ln(47r) or to the Euler constant 7 B are omitted, as they have no effect on the MS anomalous 
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dimensions. The Riemann zeta function value £(3) — 1.20206 is denoted by (3. The number 
of active flavours in the effective theory and the sum of their charges are denoted by / and 
Q, respectively. For Zyi and Z27, the charge Q u of the charm quark is retained arbitrary. 
The incoming quark charge is set to — \ in all the expressions for Z i7 , which means that they 
correspond to the down- type quark gg'7 coupling after substituting Q u = |. For the external 
up-type quark case, one should multiply all the terms in that come with no charge factor 
by —2, and then replace Q u by Qd = — §■ The same rule applies to the matrices that are 
given in Eqs. (I5.9p - fl5.lip at the end of this section. We find 
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39955153 f , 708617 f2 , 1424 A _ ( 17755 
59049 ■> 19683 ■/ ~ l ~ 6561 ^ \ 243 



80 f _ 1 oT^A 1 «i / 1 [ 925216 _ 19720 f 1 800 
81 J iz Vy "I" e \ e I 656i 729 i "r 729 

t 8 . f 2 _ f 1468 , 28 f \ 77I _i_ «1 f J_ [ 82789640 
87^ \ 27 ^~ 9 J J ^ } ^ e \e 2 [ 59049 

i f _|_ 8 f 2^ ?7l 1 1 [ 20027615 _ 39955153 f 1 701 
J ' J J t - [ 6561 59049 J i- 1S 



154280 /• 1 448 f'. 
6561 J 2187-' 

( 37826 _ 1885 f , 
V 27 9 J ' , 

28 ^2\ Tjl _ 2247378701 _ 766276 a , / 195301451 _ 21056 
) J )^\ 531441 6561 S3 ~r \^ 354294 2187 



1432 f3 , 
6561 J ^ 



(2- 

V 27 , 

" 20027615 _ . 

6561 59049 " " DODl" \ 

_ 21056 a \ f _ f 125173. 
2187 S3^| J ^ 59049 

2 a ^ £ 74 f2] 7=)\ 
81^ 



7375087 
3888 



40664 a 
^7~S3 



1251734 
590" 

q}, 



a 2 



( 62269 , 142 a \ f 

\ 324 ' 27 S3^| J 

( 1906 _ 55 A 1 ^ [1 / 873038 _ 9530 r , 695 r 2\ _ 726064 _ 410713 f _ 6031 f2 
V 81 27-/ ,j e [e \ 2187 243 ■> 486'' / 729 " 17496 ■> 2916-/ J 

J_ ( 131785021 1 448189 f _ 380411 r2 , 80 x3\ , 1 / 184494545 _ 940982911 /• , 3173417 x2 
e 2 V 19683 ' 486 ^ 8748 ■> ' 81 ^ / e V 8748 629856 ^ 104976 ^ 

2821 ^3\ _ 97818607289 _ 4146212 a , / 3023834561 , 309493 a \ £ _ ( 3278851 _ 175 A \ f2 
972 J J 5668704 2187 S3 "r ^ 3779136 ' 1458 S3yj ^ 787 32 324S3yj 



_ 997 K 
4374^ 



(5.4) 



The n-loop l/e" poles vanish in the above renormalization constants, which is a consequence 
of the lack of one-loop contributions to Z i7 and Z i8 . 

The operators (15.21) differ by a multiplicative factor of 1/5 from the bare versions of Q-j 
and Q§. We have used X; in the actual renormalization constant calculation. Each n-loop term 
is then of order a n , which makes the simple matrix relation between renormalization constants 
and anomalous dimensions applicable. Such a relation up to three-loops has been given in 
Appendix B of Ref. [12]. Extending that result to four loops in the same notation (but with 
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.Oh 



7 



0) yields 

;ll-02 



+ 
+ 
+ 



12 n-01-11 
01-11-01 



2a n a 01 -Ak ll cl 



11-01 



+ 9 C 



6 a 



13 



2a 1L a™ + 2a ul a ii a Ui + 2a n (a 01 ) 2 + 2(a°W 1 + 4« n (a 01 ) 2 - °- n " 02 c - 12 " 01 



4a 12 a 01 



2a 01 a 12 - 4a 02 a n 



5/? a 



5 



id 14 - 6a 13 a 01 



2a 01 a 13 - 6a 03 a n - 2a n a 03 - 4a 12 a 02 - 4a 02 a 12 + 4a 1 



02 a 01 a n 



4a 02 a n a 01 +2a 01 a 12 a 01 + 2a n a 02 a 01 + 2a 01 a 02 a n 



2a 11 a 01 a 02 + 2a 01 a n a 02 + 4a 12 (a 01 ) 2 



2(a°W 2 



-2a 01 a n (a 01 ) 2 



-2fa 01 ) 2 a 11 a 01 -2a 11 fa 01 ) 3 



•2fa 01 ) 3 a n 



16a 02 / t 12 + 8a 02 a 01 K 11 +4a 01 a 02 K 11 



(a 01 ) 2 K 12 



01^3 K U 



+ 0(g 



12a 03 K n 

10\ 



12a°V 13 



(5.5) 



After calculating the ADM 7 from Eq. (15.5p . we multiply the dipole operators by 5, which 
brings us to the normalization of Qr and Qs as in Eq. ( 12. 2ft . Next, we pass from 7 to 7 eff , which 
amounts to performing a simple linear transformations that stems from Eq. (12. 3p . 

The obtained complete expressions for ^ n "> cS {n = 0, 1, 2) in the b — >• 57 case (/ = 5, Q = |, 
Q« = §) are as follows: 
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486 
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172471 
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251680 
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81 


23836 
81 


6106 
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729 
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3296257 




243 




243 


243 
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2192 
81 
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(5.6) 



(5.7) 



See Eq. (38) of Ref. [7]. 
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2821 
9 
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81 


9274 
81 
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27 
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324416 
27 
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27 


31420 
9 
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81 


166432 
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81 
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27 
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For arbitrary values of /, Q and Q u , the blocks [n = 0, 1, 2) read 



243 3^ u 
— H + SQ V 



176 
81 

_88 16 

243 81 
6272 
81 



173 
162 

70 
27 
14 
27 



./ 



L 243 



~^f + 48Q 



81 

12614 _ 
2187 
2332 



74 

81 
1736 

27 
!37i 

81 



49 



54 f 

+ 36/ 

+ 160 f 
T 27 ■> 



(5.8) 



(5.9) 



_6^f _ 374(0 , 2_fQ 
2187 ^ 27 ^ u ^27-> ^ u 



5832 + 5832 J 



f + q Qu gfQ-b 



729 1 729^ 1 9 ^ u 9 
97876 4352 j _ 112 



10577 
486 
42524 



917 
972- 



729 729 
70376 _ 15788 f , _32_ e 2 _ 1407) 
2187 2187 J ' 729 J 9 ^ 



_ 2398 
243 243 



1764752 
729 



65408 f _ 313677 
729 J 3 



4193840 
L 2187 



2187 



\ f I 896 j 
J 729^ 



1136 
9 



Q 



ffQ 



159718 
729 
2281576 
243 

3031517 
729 



+ 



./ 
/ 

39719 f 253 £2 

5832 J 486^ 

15431 f _ 6031 fi 

1458 J 486 J 



, (5.10) 
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11045 fn I 316 f 2r\ , 3695 7T 
729 J V" ' 729 J ^ u ' 486 ^ 

_ 369577 
81 



77506102 _ 875374 r , 560 x2 _ 9731 ,0 , 1104! 
531441 177147 ^ 19683 ^ 162 ^ M ' 729 

15463055 I 242204 r 1120 f 2 i 55748 f) 33970 J,Q 632 f2,Q 
177147 " r 59049 ^ 6561 J ' 27 ^ M 243 ^ ^« 243 ^ ^ u 

102439553 _ 12273398 f , 5824 r2 , 26639 77 §_ fTT 

177147 59049 ^ 6561 ■> ^ 81 ^ 27^^ 

2493414077 _ 9901031 r , 243872 p _ 1184 p _ 49993 7T , 305 xTJ 
1062882 354294 J ' 59049 ■> 6561 ■> 972 ^ 27 ^ ^ 

8808397748 _ 174839456 £ , 1600 f2 _ 669694 77 , 10672 rjj 
177147 59049 ■> ' 729 ■> 81 27 ^ ^ 

7684242746 _ 351775414 /• _ 479776 r2 _ 11456 x3 i 3950201 77 
531441 177147 J 59049 ^ 6561 ^ ' 243 ^ 



421272953 
1417176 

98548513 _ 
472392 

3205172129 _ 
472392 

6678822461 
2834352 

29013624461 
118098 

72810260309 
708588 



_ 8210077 f _ 1955 f2 
472392 J 6561 ■> 

5615165 f _ 2489 f2 

78732 J 2187-/ 
_ 108963529 e , 58903 x2 
314928 ^ 4374 J 
, 127999025 f , 1699073 r-2 , 505 f 3 
1889568 J *~ 157464 ■> ^ 4374^ 

_ 64260772 f _ 230962 f 2 _ 148 x3 
19683 J 243 ^ 27 J 

, 2545824851 r _ 33778271 f2 _ 3988 f3 
472392 J 78732 J 2187 
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112216 , 728 f , 25508 r\ 
6561 729^ ' 81 

365696 _ 1168 e _ 51232 ,p 
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3508864 
2187 



64 f r\ _ 100 77 
81 J V« 27 ^ 
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27 J • '• 9 

1904 p _ 1984 
243 J 9 
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123543040 _ 207712 f , 128 j2 24880 77 _ 640 ^77 
2187 243 ■> 27 J 9 ^ 9 J ^ 

7699264 , 2854976 f _ 12320 z-2 _ 108584 77 _ 1136 fjj 
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20 x2 
9 



953042 _ 10381 f 

2187 486 
607103 _ 1679 f 
729 81 ^ 

1597588 1 13028 x 

729 t 81 J 
2312684 1 128347 r , 920 f 2 
2187 729 J 81 J 

69359224 _ 885356 f _ 5080 f2 
729 81 J 9 J 

61384768 _ 685472 f , 350 f 2 
2187 729 ^ 81 J 



(5.11) 



The /- and charge-dependence of the blocks for n — 0, 1,2 can be found in Eqs. (9)- 
(11) of Ref. [10]. The blocks are charge-independent, while their /-dependence is given in 
Eqs. (36)-(38) of Ref. [9]. 



6 Numerical Analysis 



Let £> NNLO denote the NNLO value of the B — > X s 7 branching ratio, and the superscript 
"4L— >0" indicate that the four-loop matrix is set to zero in the calculation of a given 
quantity. Normalizing the difference between £> NNLO and #nnlo ^° ^ e result, one finds 



B 



NNLO 



»4L->0 
NNLO 



v (2)eflf. 



5 



LO 



7T 



8C^ cff (^) 



(6.1) 
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4 


5 


6 


7 


8 




14 

23 


16 

23 


6 

23 


12 

23 


0.4086 


-0.4230 


-0.8994 


0.1456 


hf 


2.2996 


-1.0880 


-0.4286 


-0.0714 


-0.6494 


-0.0380 


-0.0185 


-0.0057 




0.9075 


-0.7550 


0.6971 


0.0166 


-0.7345 


0.1784 


-0.2608 


-0.0493 




-13.9117 


11.4063 


3.4986 


-0.0619 


-2.4429 


2.1210 


-0.7558 


0.1464 



Table 1: The numbers a« and h\ n ^ that occur in Eqs. (16.61) and (16. 8p . 

where fib is the low-energy scale at which the RGEs evolution is terminated, while Cj n ^ eS are 
defined by the perturbative expansion 

^) = E^)^ n)C V)- (6-2) 

n>0 

To evaluate the r.h.s. of Eq. (16. ip . we need to know the non- vanishing LO initial conditions 
for the Wilson coefficients at the electroweak scale fiQ. They read [19] 

Cf M = 1, (6.3) 

_ f0 ), \ 3x 3 -2x 2 -8x 3 -5x 2 + 7x tn A . 

= 4^lF lnX+ 24(s-l)3 ' (6 ' 4) 
(o), v ~3x 2 -x 3 + 5x 2 + 2x 

° 8 (/io) = + 8(x-l)3 ' (6 - 5) 

where x = m 2 (//o)/M^. Substituting them into the analytical solution to the RGEs (see, e.g., 
Section 3.3 of Ref. [18]), one obtains 



t(2)efif. 



/l (2) r? a 1 +2 + h W v a 2 +2 + ^8^ ^(2)^ 



8Cf )cff (/i b ) r^cfVo) + § (t^ - tfVo) + ELi *4V ' (6 ' 6) 

where 77 = a s (/io)/«s(A f b)- F° r the numerical evaluation of 77, we always use the four-loop 
RGE for (x B (ji). The numbers cij and ^■ n ' ) are given in Table [H They are found from the 

diagonalization of ■y^ cS (15.61) as well as the elements of (15.71) and 7( 2 ) cff (15.81) . Although 

^(o)eff can be 

diagonalized analytically, the decimal approximations presented in Table [1] are 
much more convenient to use. Note that X^=i h\ n ^ = 0, which follows from the initial condition 
for C7 at 77 = 1 and from the off-diagonal position of in ■y^ eS . 

The r.h.s. of Eq. (16.61) is shown in the left plot of Fig.[2]as a function of rj for x = (162/80. 4) 2 . 
The middle plot in the same figure presents this quantity as a function of fib for fi = 160 GeV 
and for the remaining input parameters as listed in Appendix A of Ref. [4] . The //^-dependence 
of the complete correction (16. ip is shown for the same parameters in the right plot of Fig. [2j 

It is interesting to compare the four-loop effect to the analogous three-loop one at the Next- 
to-Leading Order (NLO). Let the superscript "3L— s>0" indicate neglecting the three- loop matrix 
Then 



t(l)eff/ 



(6.7) 



12 



and 



(o), 



(6.8) 



The ratios (16.71) and (I6.8P are shown in Fig. [3] as functions of r\ and 





Hb [GeV] 



A*6 [GeV] 



Figure 2: The r.h.s. of Eq. (16.61) as a function of rj (left plot) and fi b (middle plot). The right plot 
shows the relative NNLO correction (16.11) to the branching ratio as a function of fif,. 




0.2 0.4 0.6 0.8 




fib [GeV] 



fib [GeV] 



Figure 3: The r.h.s. of Eq. (I6.8j ) as a function of r\ (left plot) and (middle plot). The right plot 
shows the relative NLO correction (16. 71 ) to the branching ratio as a function of 



Comparing the right plots in Figs. [2] and [3l one can see that the two corrections are similar 
in magnitude. The suppression of the NNLO one by an additional factor of a s (/ib)/7r tends to 
get compensated by larger values of h\ as compared to . However, as the remaining plots 
show, the ratios f ]6.6f) and ( 16.81) do not contradict the naive expectation of being "numbers of 
order unity" . Thus, there seems to be no deep reason behind the observed relative smallness of 
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the NLO effect and the relatively large value of the NNLO one. It should be remembered that 
the other NLO corrections exceed ~ 30% while the other NNLO ones stay within ~ 10% [3]. 

It is evident from Fig. [2] that the /independence of the NNLO correction ( 16.11) for fi b G 
[1, 10] GeV originates almost entirely from the overall factor a'j.(fib)- The ±3% higher-order 
uncertainty in £> NNLO that was estimated in Refs. [3,4] takes this /^-dependence into account. 
The central value of £> NNLO was calculated there for \i h = 2.5 GeV. At that scale, our present 
correction amounts to around —4.4%. It is worth noting that ^nlo calculated along the 
formulae of Ref. [4] turns out to be much less /^-dependent than the more complete £> nnlo . 

The actual numerical calculation in Refs. [3,4] included contributions from B$ but not from 

(2) 

the very recently found B> s . This fact is practically irrelevant for the above discussion. The 

(2) (2) 

numerical effect of B\% on £> NNLO is about 10 times smaller than the one of B\ 7 , and has the 

(2) 

same sign. The small contribution of B\ 8 will be taken into account together with other similar 
ones in a future upgrade of the phenomeno logical NNLO analysis. 



7 Summary 

We have evaluated the complete four-loop anomalous dimension matrix B^ that is necessary 
for determining the effective FCNC couplings qq lr f and qq'g at the NNLO in QCD. The 
results are presented in a form that applies to any external flavour case, with an arbitrary 
number of non-decoupled quarks. Adding an essential contribution to the NNLO QCD analysis 
B — > A s 7 has been the main purpose of our calculation. The obtained 0(a 2 s (/x&)) correction 
to the branching ratio of this decay amounts to around —2.9% for /i b = 5 GeV, and —4.4% for 
Hb = 2.5 GeV. 
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Appendix A 

This appendix contains solutions to the RGEs for all the Wilson coefficients that matter for 
b — > 57 at the NNLO before including higher-order electroweak corrections. A generic solution 
to the RGEs flZ3| reads 

C cS (fx b ) = U(fi b ,fi )C cS (fi ). (A.l) 
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k 


I 


n = 


n = 1 


n = 2 


1,2 


1,2 


2 


none 


3,4,5,6 


1,2,4 


2 


none 


7 


all 


1,2,4,7,8 


2,7,8 


8 


1,2,4,8 


2,8 


none 



Table 2: Indices k and Z of the relevant for n = 0, 1, 2. 



The block-triangular structure of the matrix U corresponds to that of (^y cS ^j , taking into 

account that the blocks and in Eq. (12.61) have a block-triangular form, too. In 
particular, U&? = 0, and Um = for k < 2 and / > 2. 

Using the perturbative expansions for Cf (fi) (16.21) and 

f>(^, / io) = E« ri (/ u o)f> W , (A.2) 

n>0 

one easily finds that 

= r ? [^°)(7W cfr ( /Uo ) + ^ 1 )C7( ) cff ( / io) 



<? (1)cff M 



rj 



(A.3) 
(A.4) 

(A.5) 



The matrices are functions of r] only 



u ki 



(A.6) 



3=0 i=l 

The powers Oj have been given in Table [TJ Only a relatively small set of non- vanishing U kl 
S7 at the NNLO, so long as 

= 0, for Z = 1,3,4,5,6, 

= 0, for I = 3, 5, 6, 



(n) 



is 



relevant for 6 

?(0)eff 
t(l)eff 



^ U;en (/^o) 



(A.7) 
(A.8) 
0, 1, 2 are 

collected in Table [2j The "magic numbers" that occur in the expressions flA.6[) for these 



which is the case in the SM. The corresponding values of the indices k and Z for n 

l kli 



uf^' are given in Tables [3HSJ 

As far as ([1q) are concerned, their SM values for % = 1, . . . , 6 can be found in Section 2 

of Ref. [16]. The coefficients cf from that paper combine with a relative minus sign to our 



d n \ i.e. a 



(n) 



i 

For i 



of Ref. [17]. They combine to our according to the relation C- n) = C\ 
upper index gets shifted because Qi and Qg in Ref. [17] were normalized as X-j and X$ in 
Eq. (15.21) here. Once C\ n \iio) are evaluated, Eq. (12.31) should be applied to obtain c\ n ^ eS (//o)- 



7, 8, one can find the corresponding in Section 6 



tt(n+l) 



c(n+l) 



The 
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